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A YOUNG SUPER STAR CLUSTER IN THE NUCLEAR REGION OF NGC 253 
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ABSTRACT 

We present observations of a massive star cluster near the nuclear region of the nearby starburst galaxy NGC 
253. The peak of near-infrared emission, which is spatially separated by 4" from the kinematic center of the 
galaxy, is coincident with a super star cluster whose properties we examine with low -resolution (R ^ 1,200) 
infrared CTIO spectroscopy and optical/near-infrared HST imaging. Extinction, measured from [Fe II] lines, is 
estimated at Ay = 17.7±2.6. The age of the cluster is estimated at 5.7 Myr, based on Br7 equivalent width for 
an instantaneous burst using Starburst99 modeling. However, a complex star formation history is inferred from 
the presence of both recombination emission and photospheric CO absorption. The ionizing photon flux has 
a lower limit of 7.3±2.5 x lO''-' s"', corrected for extinction. Assuming a Kroupa IMF, we estimate a cluster 
mass of ~ 1-4;|;q 5 x 10^ Mq. We observe a strong Wolf-Rayet signature at 2.06 fim and report a weak feature 
at 2. 19 fim which may be due to a massive stellar population, consistent with the derived mass and age of this 
cluster 

Subject headings: galaxies: star clusters — galaxies: individual (NGC 253) 



1. INTRODUCTION 

Starburst galaxies are an important window to earlier eras 
of star formation in the Universe. Characterized by a prodi- 
gious star formation rate occasionally in excess of ^ 100 Mq 
yr"' (jlCennicutt 1998), these galaxies are using up their avail- 
able gas at a rate that is not sustainable over a Hubble Time. 
Star forming regions in local starburst galaxies (e.g. M82 and 
NGC 253) resolve into dense "super star clusters" (SSCs), 
which are the most massive (~ 10^ ^o) example of clus- 
tered star formation. These SSCs, w ith core stellar densities 
in excess of lO'* pc"^ ( ijohnsonl2005h . represent star formation 
in an extreme environment. Estimates of cluster initial mass 
functions (IMFs), sizes, and masses are necessary to place 
SSCs in the context of Young Massive Clusters (YMCs) ob- 
served in both the MiUcy W a y and the Lar ge Magellanic Cloud 
(iBrandl et al.liT999ll 19961: iWalborn et al.ll2002l) . SSCs afford 
the opportunity of studying a large sample of coeval stars, in 
an environment analogous to that of the high-redshift Uni- 
verse. 

SSCs have been observed in a variety of star-forming 
galaxies, including NGC 4038-9 ["T he Antennae"] 
dWhitmore & Schweizeill l995'), NGC 12 75 (iHoltzman et all 
119921). NGC 1569 ([O 'Connell etaTI 1 1 994. and M82 
jMcCradv & Grahaml2007ii:iO'Connell et al.ll I995V SSCs are 
compact objects with radii ~ I - 5 pc, well fit by both King 
(1962) and Elson, Fall, & Freeman (1987) density profile s 
(iMcCradv et al.ll2003l iHuxor et al.ll2005l iLarsen et al.ll200l . 
The King profile was originally established for globular clus- 
ters, and it has been suggested that SSCs represent analogues 
to globular cluster precursors. Developments in this field 
hinge on determining the IMF of stellar clusters; a top-heavy 
IMF, characterized by an overabunda nce of massive stars , 
has been observe d in clusters in M 82 jMcCradv et al.ll2003h 
and NGC 1705 (ISternberg| II998'). An IMF biased towards 
high-mass stars would be prohibitive in forming a globular 
cluster as the majority of stars would leave the main sequence 
before reaching the advanced age of globular clusters (^10 
Gyr). 

The proximity of local starburst galaxies enables obser- 
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vations on the size scale of SSCs^ 

(3.9±0.37 Mpc, where 1" - 19 pc:l Karachentsev et al. 20031 
example of an archetypal starburst galaxy in the Sculptor 
group (a - 00H7™, S - -25° 17'). The starburst nature 
of NGC 253 is thought to result from th e presence of a 7 
kpc b ar that funnels gas into the nucleus (E ngelbracht et all 
119981) . Compact radio sources corresponding to H II regions 
and [Fe II] sources tracing supernova remnants (SNRs) have 
been observed in the inner star burst disk (radius 15" - 20"; 
280 - 380 pc) of NGC 25 3 (lUlvestad & Antonuccil fT997l 
Alonso-Herrero et al. '2003). Based on size measurements 
of SNRs, Lenc & Tingav (2006) infer a supernova rate of > 
0.14[i//10'*] yr"' in NGC 253, where ly (- 10,000) is the as- 
sumed supernova expansion rate in km s"' . 

NGC 253 harbors multiple discrete sources which have 
been observed at several wavelengths, including approxi- 
mately 60 compact radio sources (Turner & Ho 1985). The 
radio nucleus of NGC 253 at 2c m, TH2, is thought to be the 
kinematic center of the galaxy ([Turner & HollI985l) . Due to 
its high brightness temperature (r ~ 10^ K) and low spec- 
tral index (a ~ -0.2, where Si, oc i^"), TH2 is l ikely a syn- 
chrotron source as oppos ed t o an H II reg i on dNakai et al.l 
II995h . iMohan et al.l(l2002h and lWeaver et al.l(l2002h associate 
this source with a heavily obscured AGN. 

The peak of near-infrared emission in NGC 253 is spatially 
separated from TH2 by nearl y 4" (~ 76 pc) and emits 10% of 
the galaxy's bolometric flux ( iMohan et al.l2002h . The near-IR 
peak is thought to be a super star cluster due to its strong re- 
combination line flux (IWatson et al.lI996l) and the presence of 
[Fe II] emissio n, which is associated with SNRs deriving from 
massive stars (lAlonso-Herrero et al.ll2003h . lEngelbracht et al.l 
(1998) observed the velocity-broadened CO(2,0) feature of 
this cluster and inferred a dynamical mass of 3.9±0.9 x 10^ 
Mq within a 7 "5 radius. It should be noted that the aperture 
size used by lEngelbracht et al.l (119981) is substantially larger 
than that of a typical SSC (< I"); this inferred mass is conse- 
quentl y an overestimate of the true cluster mass. IWatson et al.l 
d 1 9961) examined the cluster and inferred a core radius of 0.7 
pc, assuming a Gaussian density profile and a distance of 2.3 
Mpc. 

Here, we investigate the properties of the super star clus- 
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TABLE 1 
OSIRIS LINE FLUXES. 



Line Flux (10 erg s ' cm ~) 



Pa/3 


2.6 ± 0.2 


Br7 


7.5 ± 0.4 


[FeH] AL644 


3.4 ±0.3 


[Fell] A 1.257 


1.0 ±0.2 


He I 


2.6 ± 0.3 


H2 1-0 S(l) 


1.3 ±0.2 



ter near the nucleus of NGC 253. Specifically, we examine 
this object's extinction, ionizing photon flux, colors, age, and 
mass using both high-resolution imaging and spectroscopic 
data. In ^ we introduce the data, ^presents a calculation 
of the extinction, as inferred from [Fe II] lines. Cluster proper- 
ties and stellar population modeling are presented in ^ Dis- 
cussion and conclusions follow in ^and ^ respectively. 

2. OBSERVATIONS AND DATA REDUCTION 

We present two complimentary data sets: low-resolution 
near-infrared spectroscopy from the OSIRIS spectrograph on 
the 4m Cerro Tololo Inter-American Observatory (CTIO) 
Blanco telescope in Cerro Tololo, Chile and optical & near- 
infrared imaging from the Hubble Space Telescope (HST). 

2.1. OSIRIS Spectroscopy 

Spectroscopic data were obtained with the OSIRIS spec- 
trograph on the 4m CTIO Blanco telescope at Cerro Tololo, 
Chile on 2002 July 25-26 by Michael Liu. OSIRIS operates 
between 0.9 - 2.4 ^m with a 1024 x 1024 HgCdTe array with 
18.5 nm pixels. The spectra were taken with a 1 "2 x 30" slit, 
with a cross-dispersing element so that the J, H, and K bands 
were covered simultaneously in adjacent orders. 

Data were obtained at three separate spatial positions along 
the nuclear region of NGC 253, with the first observation at 
the near-IR peak and subsequent pointings moved to the north 
by approximately one slit width. In all cases, the slit was ori- 
ented East-West and on-source integration time was roughly 
20 minutes. The data examined in this paper are of the near-IR 
peak. 

The spectra were rectified using the IDL routine Redspec' 
and dark subtraction, flat fielding, cosmic -ray rejection, and 
optimal extraction were performed in the standard manner. 
Wavelength calibration was achieved with helium, argon and 
neon gas-discharge lamps internal to OSIRIS and telluric cor- 
rections were completed using the AlV standard HD 12206, 
with spline fits to remove the hydrogen absorption lines. We 
measured the spectral resolution of the CTIO data as R = A 
/A A = 1,200 or 250 km s~' (derived from the width of un- 
resolved arc lines). Flux calibration was achieved using the 
spectrum of HD 12206. Integrated line fluxes and equivalent 
widths were calculated by subtracting off a best-fit linear con- 
tinuum and directly summing the line profile. 

2.2. HST Imaging & Photometry 

The typical size of SSCs is '--^ 1 - 5 pc, which corresponds 
to 50- 250milliarcseconds (mas) at the distance of NGC 253. 
Measurement of the physical properties of SSCs therefore re- 
quires very high angular resolution. Furthermore, SSCs are 

' http://www2.keck .hawaii.edu/inst/nirspec/redspec/index.html 



often heavily shrouded by dust. Infrared measurements are 
therefore preferable to optical, as Ak ^ O.lAy. 

We made a mosaic of the central 770 pc of NGC 253 from 
dithered HST Near Infrared Camera and Multi-Object Spec- 
trometer (NICMOS) archival images taken with the F160W 
and F222M filters, in addition to an archival Advanced Cam- 
era for Surveys (ACS) Wide Field Channel image taken with 
the F8 14W filter^. Several SSCs and a prominent dust lane are 
visible in the mosaic (Figure[T]). We also examined the super 
star cluster at the near-IR peak in five passbands^: F160W, 
F222M, F187N, F190N, and Paa, and aperture photometry 
was performed in three NICMOS passbands, where the cam- 
era used is indicated in brackets: FllOM [NICl], F160W 
[NIC2], and F222M [NIC2]. Due to the variable sky back- 
ground in the vicinity of the cluster, standard photometry 
packages were avoided. We chose to develop in-house IDL 
routines to perform the photometry; we outline the procedure 
here. 

In order to remove background emission and isolate the 
cluster light, we placed a 2 pixel wide annulus at a distance of 
10 pixels from the cluster center and fitted a plane to the points 
enclosed in the annulus. Using photometric conversions from 
the HST Data Handbook"*, we summed the light within a ra- 
dius of 11.5 (6.5) pixels for the NICl (NIC2) observations 
and applied an aperture correction (AC) to account for resid- 
ual cluster light not enclosed by the adopted radius, where the 
AC is given in parentheses following the filter name: Fl lOM 
(1.13), F160W (1.19), and F222M (1.22). We converted to 
physical flux units using Vega as a zeropoint and obtained the 
following photometric results: mi6o-'W222 = 1.4±0.13, mno- 
mi6o = 2.1 ±0.25, and m222 = 11.3±0.07. Errors were esti- 
mated through Monte Carlo simulations; multiple artificial 
clusters of the same magnitude as the real cluster were ran- 
domly placed on the image and their photometry recovered. 
The standard deviation of the photometric distribution was 
adopted as the error. 

3. EXTINCTION IN THE NUCLEAR REGION 

Variable extinction in NGC 253 is evident by the large, 
clumpy dust lane cutting through the nuclear region from 
the northeast to the southwest. [Fe II] lines populate the 
near-infrared spectrum (Figure|2l) and are consequently useful 
probes of extinction. The features at 1.257 /zm and 1.644 /im 
share the same upper atomic l evel and have an intrinsic lin e 
ratio of A1.257/A1.644 = 1.35 (iNussbaumer & Storevl[T988h . 
The ratio A\IAv is calculated given equations 2a and 2b in 
ICardeUi et all (Il989l) . assuming = 3 . 1 : 

Al26^^^278; ^1^=0.182 (1) 

Ay Ay 

In order to calculate the extinction to the nuclear region of 
NGC 253, we employ the general expression relating flux and 
magnitudes [m\-m2 = -2.51og(fi/f2)] and re-arrange to solve 
for the flux ratio: 

^'■26 Mm _ -^Q-().4A,.2f, ^'-64 _ -^q-0.4A,.64 ^J) 

[Fl26 ^lm\ Q [Fl64 ^lm\ g 

^ These data are from the Multimission Archive at STScI (MAST) 
database. The original data, from Pis M. Rieke [NICMOS] and J. Dalcan- 
ton [ACS], were obtained in August 1998 and September 2006, respectively. 
The exposure times were 191, 511, and 1534 seconds, respectively. 

^ These are fully calibrated and reduced NICMOS data obtained from 
MAST, also from PI M. Rieke's August 1998 program. 

http://www.stsci.edu/hst/HST_overview/documents/datahandbook/ 



Kornei & McCrady 



3 




00''47"'34.5' 34.0' 



33.5' 33.0' 
a (J2000) 



32.5' 



32.0' 



Fig. 1 . — Logarithmic-scaled mosaic of the nuclear region of NGC 253, with imaging in the K band from HST NICMOS (red), H band from HST NICMOS 
(green), and / band from HST ACS (blue). The white outline indicates the extent of the image in the H and K bands, the square represents the near-IR peak (the 
region of the observations), and the cross indicates the position of TH2, the brightest source at 2cm and the assumed kinematic center of NGC 253 (Turner & 
119851) . Note the prominent dust lane and point-like super star clusters extending diagonally in the image from the northeast to the southwest. [See the electronic 
edition of the Journal for a color version of this figure.} 



These equations are then divided by one another to yield an 
expression which includes the observable on the left hand side 
(0.28±0.07) and the fiducial ratio (1.35) as the first term on 
the right hand side: 



Fl.26 
Fl.64 /jm 



Fi, 



26 /iin 



Fl.64 /j.m 



10 



-0.4/4 1 



10 



-0.4Ai. 



(3) 



The above expression is solved in terms of A\26 fim and 
^1.64 /jm, which are in turn related to Ay by equation [T] We 
calculate Av = 17.7±2.6. 

While an estimate of the extinction enables corrections 
to be applied to starburst parameters ( ^4.2b . it is difficult 
to constrain Ay globally due to variability along differ- 
ent li nes of sight from t he uneven dust distribution (Fig- 
ure |3j IWatson et al.lll996h . Literature values for the extinc- 
tion to the nuclear region of NGC 253 r ange from Ay = 2 
(iMoorwood & Oliva 1988)^ to Ay = 24 (ISams etal.iri994 . 
The discrepancy between these values i s attributed both to the 
highly variable extinction in the region (lAlonso-Herrero et alj 
12003b and the different aperture sizes used by various groups. 

4. CLUSTER PROPERTIES 

In order to investigate the properties of the super star 
cluster at the near-IR peak of NGC 253, we used the 
Starb u rst99 population synthesis m odel suite (iLeitherer et alJ 
Il999t IVazquez & LeithereJ l2005h to create a hypothetical 



^ The IMoorwood & Olival jl988l) value for Ay appears to be incorrect, as 
a derivation using their input parameters from their Table 3 would yield Ay 



Starburst cluster and model its properties as a function of time. 
We assumed solar metalli city, as p er studies of NGC 253 by 
ICarral et alj ([1994) and P tak et all (119971) . and modeled both 
an instantaneous burst (e.g. a "single stellar population") and 
continuous star formation. Both of these models may be over- 
simplifications of the super star cluster at the near-IR peak of 
NGC 253, as discussed below. 

The stellar IMF is often modeled as a single power-law 
across all masses (e.g. the Salpeter (1955) IMF, where dN 
oc M""dM; a = 2.35), although recent work has suggested 
that the IMF may be better represented as a broken power law 
(iMeyer et al. 2000). An example of such an IMF with a char- 
acteristic "knee" in its distribution is the Kroupa (2001) IMF; 
this broken power law has the effect of reducing the number 
of low-mass stars relative to a population described by a sin- 
gle power law. We use Starburst99 to model a stellar popula- 
tion according to the Kroupa IMF. We first investigate the age 
of the cluster using Starburst99 and the Br7 equivalent width 
measurement. 

4. 1 . Equivalent Width & Age 

Br7 (2.17 fim) is a bright recombination line indicative 
of a massive, young stellar population. As high-mass stars 
evolve much more quickly than their low-mass counterparts, 
the strength of the Br7 feature will vary over time as the out- 
put of UV photons fr om massive stars drops sharply after the 
first 5 Myr (Leitherer et al.lfT999h . Therefore, the strength of 
the Br7 feature, quantified by its equivalent width (EW), is a 
probe of cluster age. As both line flux and continuum level in- 
crease with cluster mass, EW is relatively insensitive to mass 
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Fig. 2. — J, H. and K band CTIO spectra (R ~ 1,200) of the nuclear region of NGC 253. The signatures of CO rovibrational transitions in the H and K bands 
are due to photospheric absorption in red supergiants, while the recombination lines Pa/3 and Br7 are indicative of a population of OB stars. [Fe II] lines in the 
J and H bands trace supernova activity and can be used to calculate the visual extinction (Q. The strong detection (9o") of He I at 2.06 im\ and the possible 
detection of He II at 2.19 /im suggest the presence of Wolf-Rayet stars, though higher spectral resolution observations are necessary in order to conclusively 
identify such a population. 



and is consequently a powerful age diagnostic. 

From the CTIO data (Figure |2l), we measured the equiva- 
lent width of the Br7 feature to be 25.3±1.3 A. We used Star- 
burst99 to investigate the temporal evolution of Br7 equiva- 
lent width (FigurelU. The two instantaneous burst models we 
used — one with a Kroupa IMF and one with a top-heavy 
Salpeter IMF — produced indistinguishable distributions of 
Br7 equivalent width with time (as would be expected, since 
these models are nearly identical in their distributions of high- 
mass stars). The continuous star formation model sustained a 
higher Br7 equivalent width over time than the burst models 
did. Continuous star formation may be an unrealistic assump- 
tion, as the winds from massive stars and their ultimate explo- 
sions as supernovae would be effective in clearing out natal 
gas and therefore impeding subsequent star formation. Based 
on the assumption of a single stellar population, the age of the 
cluster from Br7 equivalent width measurements is 5.7 Myr 
(main sequence turn-off mass of ^ 20 Mq). The Br7 line is 
well-resolved in the CTIO data and therefore the uncertainty 
on the cluster age is dominated by the choice of star forma- 
tion history model. Assuming a metallicity other than solar 
minimally affects the best-fit cluster age; adopting Z = 2Zq 
(O.5Z0) produces a best-fit age of 5.4 (6.3) Myr. 

Next, we investigate another defining parameter of the clus- 
ter - its mass. To do so, we first calculate the cluster's ionizing 
photon flux. 

4.2. Ionizing Photon Flux 

Measurements of SSC recombination line luminosity can 
be used to calculate the flux of ionizing photons and conse- 
quently the number of stars producing the UV photons (as- 



suming Case B recombination). The flux of ionizir ig photons 
is calculated from the Br7 luminosity, as given in iKennicuti 
(1998): 

8.2 X 10-4" Lfir^ [erg s"'] = 1.08 x 10"" Q(H°) [s"'] (4) 

where Q(H*') is the number of hydrogen ionizations per sec- 
ond and Lbi-/ is defined as the flux in the Br7 line integrated 
over the surface area of a sphere whose radius is the dis- 
tance to the source. We calculate Q(H°) = 1.0±0.2 x 10^^ 
s"', uncorrected for extinction. Correcting for attenuation in- 
creases the value of Q(H") by a factor of lO''-^'*^ " where 
Aim ^J.m = 0.12Av using the Cardelli (1989) extinction law 
(Ry = 3.1). Based on the extinction to NGC 253 inferred from 
[Fe II] lines in §4, Q(H'')con- = 7.3±2.5 x 10^^ s"'. The value 
of Q(H°) which w e derive is comparable to the result from 
lPuxlevetal.l(ll997h using radio recombination lines. Discrep- 
ancies between literature values likely arise primarily from 
the use of different distances t o NG C 253; some authors, 
includ ing I Alonso-Herrero et alj (l2003h and Engelbracht et ajQ 
(119981) . adopt a value of 2.5±0.37 Mpc (f de Vaucouleurs. 19781) 
whereas this work uses the most current v alue for the distance 
to NGC 253 (3.9±0.37 Mpc) obtained bv lKarachentsev et alj 
(12003). If we scale Q(H")corr to a distance of 2.5 Mpc, we ob- 
tain a value of 3.0±1.5 x 10^^ s"'. A secondary effect which 
may further explain differences in Q(H") is the use of differ- 
ent extinction laws. The Galactic extinction along the line of 
sight^ to NGC 253 is Ak = 0.007 and therefore the lower limit 
on Q{lf), assuming a distance to NGC 253 of 3.9 Mpc, is 
1.0±0.2 X 10^^ s"'. Table |2] summarizes the ionizing pho- 

* http://nedwww.ipac.caltech.edu/forms/calculator.html 
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Fig. 3.— Comparison of the nuclear region super star cluster in NGC 253 observed by HST NICMOS [NIC2; Of'075 pixer'] in five passbands (F160W, 
F222M, F187N, F190N, and Paa) and HST ACS [WFC; 0"050 pixer'] in the F814W passband. Each linearly-scaled image is 3"75 on a side and the image 
y axis is oriented 68° east of north. Note the Airy ring in the F222M image; the source is only marginally resolved at this wavelength and therefore reflects the 
point spread function. The higher emission on the left side of the cluster relative to the right side (especially visible in the F187N and Paa images) is in agreement 
with the placement of the large dust lane in the nuclear region of NGC 253 (Figure[T). The Paa image also shows clear signs of asymmetry, possibly indicative 
of an outflow. 
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Fig. 4. — Starburst99 model of Br7 equivalent width (EW) versus time, 
for three combinations of star formation histories and IMFs. The instanta- 
neous burst models (solid and dot-dashed lines, where the latter is coincident 
with the former) have a characteristic fall-off in Br7 equivalent width around 
the main sequence turn-off age of O stars (~ 5 Myr) while the continuous 
star formation model (dotted line) sustains a relatively constant Br7 equiv- 
alent width until an age of approximately 50 Myr The observed value of 
Br7 EW (25.3±1.3 A) has been plotted as a solid horizontal Une. Assuming 
that the SSC in the nucle ar re gion of NGC 253 experienced an instantaneous 
burst of star formation ( j|4. IK the age of the cluster is estimated at 5.7 Myr. 
Uncertainty in this estimation is dominated by model choice, as opposed to 
measurement of the Br7 equivalent width. 

ton flux derived here, assuming a distance of 3.9 Mpc, and 
compares it to that of other authors. 

Ionizing photon flux, as a proxy for the number of massive 
stars in the cluster, can be used to infer a total cluster mass. 
We address this calculation in the following paragraph. 

4.3. Cluster Mass 



SSCs have been observed with masses similar to those 
of globular clusters, ranging from 10**-^ - 10^-^ Mq (e.g., 
iMassev&Hunteill 19981: iGilbert & Grahamll2007h . Here, we 



TABLE 2 
Flux of ionizing photons. 



Log Q(H") [s 



Method 



Reference 



53.86 ± 0.40 

53.57 

53.15 

53.11 

53.00 

52.93 

52.88 



Br7 (extinction connected) 

RRLs" 

Pao 

[Ne II] 

Br7 

RRLs" 

6cm thermal flux 



This Work 
Puxlev et al. (1997) 
Alonso-Herrero et al. (2003) 
Engelbracht et al. (1998) 
Engelbracht et al. (1998) 
Mohan et al. (2002) 
Turner & Ho (1983) 



" Radio recombination lines (n > 50). 

investigate the mass of the super star cluster near the nucleus 
of NGC 253 using a novel technique which relies on the ob- 
served total cluster ionizing flux. 

We integrate an assumed IMF (Kroupa) over mass, using 
the observed ionizing photon flux, 7.3±2.5 x lO" s"' (EH, 
as a constraint. This procedure requires a priori knowledge 
of th e ionizing flux of indi vidual stars - we use the values 
from [Sternberg et al] (l2003h . We included only stars of mass 
20-90 Mq in our analysis, as stars below 20 Mq contribute a 
comparatively insignificant ionizing flux. 

With known stellar ionizing fluxes from [Sternberg et alj 
(|2003) in hand, we made an arbitrary, first-guess estimate of 
the cluster mass (10^ Mq). Next, we calculated the distri- 
bution of stellar masses necessary to produce this total mass, 
assuming a Kroupa IMF. The number of stars in each mass bin 
were then multiplied by their respective ionizing fluxes and a 
total ionizing flux was extracted. This value was compared 
to the observed ionizing flux and the total mass of the clus- 
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Fig. 5. — Near-infrared color-color plot from Starburst99 modeling. The 
tracks assume a Kroupa IMF, and an instantaneous and continuous star forma- 
tion history, respectively. The photometric output of Starburst99, in Bruzual 
filters, has been transformed to the NICMOS photometric system using the 
method of Marleau et al. (2000). Cluster ages, in Myr, are indicated by black 
squares. The measured photometry of the cluster, with Icr errors, is indicated 
by the cross, and a reddening vector assuming a Cardelli (1989) extinction 
law is extended. Through this technique, an extinction of Ay ~ 12 is in- 
ferred. Given the differential extinction towards the cluster, the uncertainty 
on the photometry translates into a cluster age between 1-4 Myr. 

ter was iteratively adjusted until the model and observed total 
ionizing flux were in agreement. We obtain a best-fit cluster 
massof 1.4![|;^ x lO'^M©. 

4.4. Cluster Photometry 

Cluster photometry is revisited in light of Starburst99 mod- 
eling. Extinction to the region was previously determined 
(|3]) using observations of [Fell] lines; we now independently 
examine the extinction level using Starburst99 modeling and 
HST photometry (©. We made a near-infrared color-color 
diagram for two model stellar populations; both assume a 
Kroupa IMF and one is characterized by an instantaneous 
burst while the other has a continuous star formation history. 
The output photometry from Starburst99 was converted from 
the Bruzual (1981) JHK system to NICMOS magnitudes us- 
ing the method of iMarleau et al.l (2000), and the observed 
HST photometry plotted with error bars (Figure|5]l. The inter- 
section of the reddening vector (calculated using the Cardelli 
(1989) extinction law) with the color-color track implies an 
extinction of Av/ ^ 12 and a corresponding age of r^ 1 - 4 Myr 
Given the cluster's high mass and young age, the spectral sig- 
natures of a population of Wolf-Rayet stars are expected (|3. 

5. DISCUSSION 

The nuclear region of NGC 253 is heavily extinguished, 
as seen by the dust lane running from the northeast to the 
southwest in Figure [T] Infrared observations are therefore 
preferable to optical observati ons, not only for the advantage 
in attenuation {Ak ^ 0.1 lAy; ICardelh et 31111989 ), but also 
because of the wide range of nebular and photospheric lines 
which are most pronounced in the infrared. 



The super star cluster at the near-infrared peak of NGC 253 
appears to be a complex system consisting of several epochs 
of star formation. CO absorption features characteristic of 
photospheres of red supergiants (RSGs) are seen in the H and 
K band spectra (Figure |2]i. Stellar winds from these massive 
stars remove remaining nascent material from the cluster and 
near-infrared light becomes dominated by the spectral signa- 
ture of RSGs, rich in OH, CO, and H2O features. RSGs are 
evolved OB stars, so their presence in the cluste r implies a 
minimum age of - 7 Myr Ceitherer etdll 19991) . The data, 
however, also exhibit strong recombination Unes due to OB 
stars, indicating the presence of a young population still pro- 
ducing UV photons. A single burst of star formation is un- 
able to explain the presence of both RSG and OB star spec- 
tral signatures and therefore either a continuous star forma- 
tion history or one punctuated by multiple bursts is necessary 
to account for the features seen in these spectra. Other possi- 
bilities which could explain the seemingly incongruous result 
of simultaneously finding both OB and RSG stars in a clus- 
ter include the presence of binary systems (and the associated 
processes of mass transfer and interactions) and line-of-sight 
superposition of multiple clusters. Furthermore, the observed 
cluster could in actuality be a ble nd of two or more bursts 
of star formation; ISeth et"an ( l2008l) present observations of a 
nuclear cluster in NGC 4244 with two distinct stellar popu- 
lations. Explanations for this "blended" cluster include peri- 
odic bursts of star formation triggered by gas accretion onto 
the cluster and merging of multiple SSCs due to dynamical 
friction. 

We used archival HST imaging of NGC 253 (Figure [3]) to 
examine the dusty nuclear environment. The large dust lane 
clearly visible in Figure[T]is also evident in the postage-stamp 
images of Figure[3]in the form of increased extinction towards 
the left (south). Furthermore, the asymmetry observed in the 
Paa image may be indicative of an outflow from the clus- 
ter. Alonso-Herrero et al. (2003) find diffuse Paa emission 
both in and around the disk of NGC 253. An outflow might 
imply that the cluster is not bound, which would reduce the 
possibility of its development into a globular cluster. Older 
clust ers are observed to be more compact than younger clus- 
ters dMenge l et al. 2005), although the apparent evolutionary 
advantage of being compact has yet to be fully explored in 
light of the selection bias against diffuse clusters (regardless 
of age). The "infant mortalit y rate" of SSCs has been esti- 
mated to be as high as 99% ( Fall & Zhangll2001h . The pro- 
posed evolutionary progression of SSCs into globular clusters 
hinges both on SSCs remaining bound for a period of several 
Gyr and SSCs having an IMF comparable to that of globular 
clusters. 

Given a young cluster with a sufficiently well-sampled IMF, 
the spectral signatures of a population of Wolf-Rayet (WR) 
stars are expected. The nuclear cluster in NGC 253 appears 
to be both massive and young enough to harbor such a popu- 
lation. The strong emission line signatures of WR stars in the 
near-infrared include He I at 1.08 /im and He II at 1.01 /im, 
in addition to weaker features at 2.06 /im (He I) and 2.19 /im 
(He II). The 2.06 /im feature (EW - 8.8±1 A) is detected at 
the 9(7 level in the CTIO spectra, and we also report a possible 
detection of the He II feature at 2.19 /xm (Figure|2]i. Follow- 
up observations with higher spectral resolution are necessary 
to conclusively indicate the presence of WR stars. 

The short lifespan of a WR stars limits the age of a single- 
burst population to a range of 3 - 8 Myr. The equivalent 
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width of WR signatures, used in combination with the equiv- 
alent width of h ydrogen reconibinati o n features, is a powerfu l 
age diagnostic ( Schaerer et al.lll999l) . ICrowther etal] (l2006h 
used emission line ratios in the J and K bands to classify WR 
subtypes in the young Galactic cluster Westerlund 1, demon- 
strating the util i ty of t he near-infrared in studying WR stars. 
iLumsden et alj (1 19941) conducted a search for WR stars in 
three starburst galaxies and failed to find any signatures in the 
K band; the 2.06 and 2.19 /im features are more susceptible 
to dilution from the continuum of older background popula- 
tions than are the WR lines in the J band. O bservations of 
WR stars in starburst galaxies have been rare (IContinilll996t 
iMas-Hesse et al.ll20odlMassey & Hunte3ll998 '). although the 
nuclear cluster in NGC 253 is a plausible candidate for host- 
ing WR stars due to its large mass and low age. 

6. SUMMARY AND CONCLUSIONS 

We have used spectroscopy and imaging to examine the 
super star cluster coincident with the near-infrared emission 
peak in the nearby starburst galaxy NGC 253. Line ratios, 
equivalent widths, and integrated line fluxes were measured 
using low resolution (R 1 ,200) infrared CTIO spectroscopy. 
Extinction, estimated using [Fe II] lines, is Ay = 17.7±2.6. 
This value was used to correct the flux of ionizing photons; 
Q(H")corr = 7.3±2.5 X 10^3 s"'. Assuming a Ki'oupa IMF, 
this flux of ionizing photons implies a total cluster mass of 
lA-^i, X 10^ Mq. We examined the color of the NGC 253 nu- 
clear star cluster using aperture photometry on archival HST 
images: mi6o-'«222 = 1.4±0.13, mwf^-mia) = 2.1±0.25, and 
ni2i2 = 11.3±0.07. We also used Starburst99 modeling, as- 
suming a solar metallicity, an instantaneous burst of star for- 
mation, and a Kroupa IMF, to infer age from the measured 
equivalent width of Br7. The strength of the Br7 line is pro- 
portional to the flux of ionizing photons, a quantity which 
depends on the number of main-sequence OB stars in the 
cluster The measured Br7 equivalent width, 25. 3 ±1.3 A, 
corresponds to a cluster age of 5.7 Myr, assuming the Star- 



burst99 parameters stated above. While this technique gives a 
rough estimate of the cluster's age, it is important to remem- 
ber that the near-infrared CTIO spectra exhibit both recom- 
bination emission and CO photospheric absorption. These 
spectral signatures of both young and older stars imply ei- 
ther a continuous star formation history or one punctuated by 
multiple bursts; a single burst is not sufficient to explain the 
cluster's star formation history. 

Future studies of this cluster to probe its stellar velocity 
dispersion would be useful in constraining its mass. High- 
resolution spectroscopy would also be instrumental in quanti- 
fying the contribution of Wolf-Rayet stars, if any, to the over- 
all cluster light. 
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